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1.0  INTRODUCTION 


1.1  Background/Scope/Program  Objectives 

Any  biological  process  is  extremely  complex,  involving  numerous  cell  types  in  a  variety 
of  states,  as  well  as  extracellular  molecules  and  cytokines.  These  biological  processes  have 
extensive  interconnections  and  complicated  feedback  loops  that  help  to  produce  stable,  self- 
repairable  systems.  In  a  healthy,  “normal”  biological  system  the  multitude  of  processes  are  in  a 
homeostatic  state.  However,  should  an  event  occur,  such  as  an  injury  or  infection,  this  stable 
state  is  automatically  modified  to  bring  about  changes  to  address  the  problem.  Eventually,  if  the 
response  is  successful,  the  injury  is  healed  or  the  infection  is  cleared,  and  the  biological  system 
can  return  to  a  healthy,  “normal”  homeostatic  state.  However,  even  when  the  response  is 
successful,  the  healing/repair  process  is  not  always  perfect  and  often  results  in  modifications  to 
the  specific  environment,  such  as  fibrosis  and  loss  of  functionality.  In  many  areas  of  the  body 
such  imperfect  healing  simply  results  in  a  blemish  on  the  skin,  a  welt  in  a  bone,  or  a  less  flexible 
joint.  However,  in  some  locations,  such  as  the  heart  or  the  eye,  imperfect  repair  can  result  in 
significant  loss  of  capability  such  as  heart  disease  or  blindness. 

A  biological  system,  as  with  any  complex  system,  can  be  broken  down  hierarchically  so 
that  different  views  on  the  system  can  be  addressed.  For  example,  the  typical  breakdown  of  the 
human  system  includes  the  digestive  system,  the  circulatory  system,  the  skeletal  system,  the 
muscular  system,  the  respiratory  system,  the  nervous  system,  etc.  These  lines  are  somewhat 
artificial,  but  drawing  such  divisions  allows  us  to  deal  with  the  extreme  complexity  of  biological 
systems.  Within  each  system  are  specific  organs  and  processes  that  perform  the  functions 
required  of  that  particular  system.  In  the  case  of  the  nervous  system,  such  components  include 
the  brain,  the  spinal  cord,  and  the  senses.  Thus,  the  eye  is  considered  part  of  the  nervous  system 
as  it  provides  visual  input  to  the  brain.  The  eye  is  an  extremely  complex  organ,  involving 
numerous  layers  of  varying  cell  types  to  perform  its  function  of  converting  photons  into  a  signal 
that  can  be  interpreted  by  the  brain  as  a  visual  image.  In  humans  this  capability  is  highly 
evolved  and  sight  is  an  extremely  important  source  of  information  about  the  world  around  us. 

The  ultimate  goal  of  any  work  in  developing  a  model  of  the  eye  would  be  to  build  a 
complete  model  of  the  cells  and  processes  found  in  the  eye,  especially  the  retina,  as  it  is  the 
primary  facility  for  sensing  light  and  converting  them  to  sensory  signals.  Building  such  a  model 
would  require  considerable  resources  and  time  to  try  to  do  all  at  once.  Addressing  such 
complexity  would  be  best  handled  by  building  a  modular  model  that  can  address  one  component 
at  a  time,  allowing  the  entire  model  to  be  built  one  manageable  chunk  at  a  time.  Because  the 
interest  of  the  AFRL711  HPW/  RHDO  is  the  effect  of  lasers  on  the  eye,  and  in  particular  the 
thermal,  photochemical,  and  photomechanical  effects  on  the  retina,  the  area  focused  on  in  this 
project  was  the  development  of  a  computational  model  of  the  complex  cellular  interactions  that 
occur  in  the  retina  of  the  eye  when  exposed  to  a  laser,  with  particular  emphasis  on  the  retinal 
pigmented  epithelial  cells,  or  RPEs.  This  model  emphasizes  the  control  systems  and  feedback 
loops  internal  to  the  RPE  that  are  involved  in  the  cellular  response  to  such  exposure,  utilizing 
any  specific  data  available  about  how  RPEs  respond  to  various  levels  of  laser  exposure.  Of 
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course,  the  response  of  RPEs  to  such  exposure  invokes  an  entire  cascade  of  events  involving  a 
typical  biological  system’s  response  to  stress,  including  the  biochemical  reactions  that  alter  cell 
states  and  the  call  for  cells  from  the  local  area  and  blood  stream  in  response  to  such  stimulation. 
Portions  of  these  interactions  were  also  modeled  based  on  their  direct  affect  on  the  RPEs.  Thus, 
the  model  also  contains  high-level  knowledge  of  the  various  cellular  responses  to  the  inputs,  and 
the  cell-to-cell  interactions  based  on  their  cytokine  and  other  chemical  output.  Input  to  the 
model  allows  time-history  of  ambient  temperature,  low-light  exposure  levels,  and  antioxidant 
availability.  A  primary  output  of  the  model  involves  levels  of  RPE  viability  over  time,  as  well  as 
the  time  course  of  many  relevant  values  available  about  the  RPE  and  the  processes  taking  place 
in  it. 


The  modeling  approach  used  in  the  work  reported  on  here  is  a  technique  called 
BioEusion®',  which  was  developed  by  the  Principle  Investigator,  and  for  which  she  is  co¬ 
inventor  on  two  patents.  This  is  a  unique,  advanced  modeling  technique  that  can  integrate  a 
variety  of  levels  of  modeling,  from  numerical  to  symbolic,  into  a  single  model,  allowing  the 
ability  to  model  large,  complex  sets  of  cellular  and  chemical  interactions  on  a  PC -based 
computer.  Models  can  be  built  from  the  sub-cellular  level  up  through  entire  disease  processes. 
The  approach  allows  for  a  systematic  collection,  organization,  and  analysis  of  the  knowledge 
needed  to  build  the  model,  as  well  as  a  modular,  structured  approach  to  building  the 
computational  model  using  a  commercial,  off-the-shelf  product. 

Current  research  in  the  effects  of  lasers  on  the  eye  have  centered  on  developing  an 
understanding  of  how  light  and  heat  affect  the  various  cells  in  the  eye,  particularly  the  RPEs,  and 
the  level  of  exposure  required  to  produce  a  certain  short-term  effect  on  the  RPEs.  Mathematical 
models  have  been  developed  (Torres  et  ah,  1993;  Till  et  ah,  2003;  Chen  et  al,  2006;  Goldberg  et 
ah,  2007)  that  represent  these  effects.  The  BioEusion  approach  used  in  the  work  reported  on 
here  produces  a  higher-level  model  of  the  problem,  modeling  a  large  portion  of  the  RPE  cell, 
including  its  “normal”  healthy  state  and  how  its  behavior  is  modified  due  to  various  levels  of 
laser  exposure.  This  involved  modeling  various  intra-cellular  processes,  including  the  internal 
biochemical  cascades  that  are  initiated  under  light  exposure,  plus  metabolic  and  redox  processes, 
and  the  effect  of  oxidative  stress  on  the  cell’s  lipoproteins,  proteins,  and  DNA.  This  approach  to 
modeling  allows  researchers  to  integrate  the  findings  from  the  mathematical  approaches  to  laser 
damage  that  are  available,  including  thermal,  photochemical,  and  photomechanical  effects  on  the 
RPEs,  with  the  knowledge  of  the  RPEs’  behavioral  response  to  such  exposure,  as  well  as  the 
cellular  interactions  occurring  in  the  RPEs’  immediate  environment.  Use  of  the  modular, 
hierarchical  approach  to  developing  the  model  allowed  work  to  focus  on  the  exploration  of  the 
effects  on  the  RPEs’  viability  and  later  to  expand  the  model  to  include  details  of  other  cell  types 
and  processes  found  in  the  eye,  and  their  reactions,  to  such  exposure. 

1.2  Modeling  Biological  Processes 

Knowledge  about  biological  processes  is  inherently  complex,  variable,  and  uncertain. 
The  levels  of  knowledge  needed  to  understand  how  a  biological  process  is  influenced  are:  1) 
sub-cellular,  chemical  interactions  that  generate  a  local  environment,  2)  individual  cell  behaviors 


^  BioFusion  is  a  registered  trademark  of  Interieukin  Genetics,  Inc. 
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within  an  environment,  3)  sets,  or  pools,  of  a  specific  cell  type  within  an  environment,  4) 
aggregations  of  different  cell  pools  whose  chemical  products  and  behaviors  create  a  particular 
environment,  and  5)  the  environment  generated  by  the  aggregation  of  cell  pools  that  manifests 
itself  at  the  system,  organ,  or  patient  level  as  the  disease’s  symptom  complex.  This  list  contains 
a  feedback  loop,  since  the  environment  listed  in  number  one  is  aggregated  to  generate  a  global 
system  environment  in  number  four,  which  in  turn  influences  numbers  two  and  three,  and  is 
manifested  at  the  system/organ/patient  level  in  number  five  as  the  symptom  complex.  The 
complex  interaction  between  the  levels  of  biological  knowledge  makes  it  advantageous  to 
represent  and  implement  the  knowledge  independently  within  a  hierarchy  to  accurately  specify 
all  relationships  between  levels  (Fink  and  Herren,  1996). 

Thus,  in  order  to  model  a  complex  biological  process,  such  as  the  processes  that  are 
initiated  in  the  RPEs  by  exposing  the  eye  to  a  laser,  sufficient  detail  must  be  provided  so  that  the 
model  can  run  as  a  simulation.  For  example,  sufficient  knowledge  must  be  available  about 
normal  RPE  cellular  behavior,  the  specific  initial  effects  laser  exposure  has  on  the  RPE 
individually,  what  response  these  cells  have  to  those  specific  initial  effects,  what  signals  RPEs 
send  externally  that  would  affect  other  RPEs  (and  other  cell  types)  in  the  area,  what  possible 
states  the  RPEs  can  be  in  based  on  their  level  of  exposure  to  the  laser,  what  environmental 
factors  might  affect  the  cells  and  their  response  to  laser  exposure,  and  what  intra-cellular 
processes  get  initiated  as  a  result  of  the  direct  exposure  as  well  as  the  feedback  from  the 
environment  generated  by  the  stimulated/damaged  RPEs. 

Methods  of  acquiring  and  representing  such  knowledge  about  biological  processes  to 
build  a  viable  model  must  surmount  a  number  of  obstacles  and  issues.  These  are  as  follows: 

1)  The  level  of  existing  knowledge  about  a  disease/repair  process  and  the  basic  supporting 
biology  tends  to  be  quite  shallow.  Researchers  often  do  not  have  information  about  exact  causal 
relationships  and  only  have  empirical  data  on  associations  between  stimulus  and  resulting 
behaviors. 

2)  Anatomy  and  physiology  of  a  particular  system  is  complex  and  variable.  Biological 
processes  have  a  large  amount  of  variation.  Not  only  do  aspects  of  anatomy  vary  from  person  to 
person  but  the  basic  biological  processes  are  greatly  modified  by  the  background,  history,  and 
genetics  of  a  particular  patient.  The  quantity  of  knowledge  that  must  be  acquired  and  encoded 
can  be  huge,  so  efforts  must  be  made  to  simplify  and  aggregate  wherever  possible  and 
appropriate.  In  addition,  space/distance  often  plays  an  important  role  in  how  various  portions  of 
the  anatomy  interact  and  influence  each  other  in  a  particular  biological  process. 

3)  Biological  systems  are  based  on  self-generated  feedback  and  control  and  are  not  static, 
but  evolving,  adapting  systems.  Biological  systems  work  because  of  the  self-regulating  feedback 
and  are,  fundamentally,  control  systems.  Eeedback  is  the  foundation  of  any  biological  system. 
Though  current  simulation  techniques  can  handle  feedback,  the  level  of  complexity  and  the 
pervasiveness  of  the  feedback  in  biological  systems  is  beyond  what  is  commonly  modeled  in 
other  more  common  model  application  areas.  The  modeling  approach  used  in  BioEusion,  which 
generates  a  simulation  over  time,  can  handle  such  complexity. 

4)  What  is  known  about  a  particular  biological  process  is  often  uncertain  and  contradictory. 
The  body  of  knowledge  in  a  particular  biological  area  will  have  gaps,  as  well  as  contradictions. 
Sources  of  information  may  only  be  tangentially  related  or  may  be  unrepresentative  studies  that 
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provide  insight  without  confirmation.  Because  of  this,  the  research  team  must  synthesize 
knowledge  from  a  wide  variety  of  sources  to  bring  to  bear  on  a  modeling  effort. 

5)  The  knowledge  on  which  the  model  is  built  is  constantly  changing  and  evolving.  The 
biological  sciences  are  constantly  uncovering  new  information  about  disease  processes.  One  of 
the  goals  of  a  BioFusion-based  computational  model  is  to  synthesize  that  knowledge  to  express  a 
comprehensive  theory  about  the  biological  process  under  study.  Once  the  model  has  been  built 
representing  the  theory,  the  theory  can  be  tested,  refined,  and  validated.  Because  the  model  is 
modular  and  hierarchical,  as  new  information  emerges,  it  can  be  incorporated  into  the  model. 

As  a  result  of  these  issues,  not  only  is  the  actual  modeling  process  difficult  in  a  biological 
area,  but  the  knowledge  acquisition  process  is  also  highly  complex.  Thus,  although  lower  level, 
highly  specific  areas  can  be  handled  using  mathematical  techniques,  at  the  higher  levels  broader, 
less  restrictive  methods  of  representation  are  needed  that  allow  for  the  complex  feedback 
interactions  and  the  large  uncertainties  in  the  knowledge.  Modeling  techniques  that  are  forgiving 
of  the  lack  and  imprecision  of  knowledge  are  needed  to  effectively  develop  useful  models  in 
support  of  biological  research.  That  is  why  we  used  a  modeling  technique  developed  specifically 
to  model  biological  processes,  namely  BioFusion.  This  modeling  approach  was  designed 
specifically  to  address  the  issues  that  make  modeling  biological  systems  complex,  providing 
methods  for  knowledge  acquisition,  as  well  as  flexibility  in  the  implementation  process  to  allow 
for  experimentation  and  theory  generation/exploration. 

1.3  Overveiw  of  the  BioFusion  Modeling  Technique 

The  BioFusion  model  development  methodology  was  used  to  design,  build,  and  validate 
a  computer-based  model  of  the  laser  effects  on  the  RPEs  in  the  eye.  Though  the  model  was 
developed  such  that  it  works  under  “non- stimulated”  conditions  as  a  baseline,  it  also  includes  the 
ability  to  provide  certain  exposures  to  the  cellular  environment,  including  laser  exposure  level, 
ambient  temperature,  oxygen  availability,  low-light  exposure  levels,  and  the  current  existing 
state  of  the  cellular  environment,  including  availability  of  various  types  of  anti-oxidants.  The  PI 
for  the  work.  Dr.  Pamela  Fink,  developed  the  BioFusion  methodology  while  at  Interleukin 
Genetics,  Inc.  (formerly  Medical  Science  Systems,  Inc.).  The  technique  received  a  patent  on 
August  12,  1997  (patent  number  5,657,255,  “Hierarchical  Biological  Modeling  System  and 
Method”).  Several  papers  have  been  published  on  this  technique  (Fink  &  Herren  1996,  Fink  & 
Herren  1997,  and  Herren  et  al.  1998)  as  well  as  on  its  use  in  a  variety  of  applications  (Fink  et  al. 
2001,  Fink  &  Morgan  2000). 

The  BioFusion  model  development  methodology  is  a  hierarchical,  iterative  approach  to 
model  development.  It  involves  gathering  together  all  of  the  knowledge  that  can  be  found  on  a 
particular  targeted  topic,  organizing  it  into  a  coherent  “theory”  of  the  process  in  question,  and 
implementing  this  “theory”  into  a  computer-based  model.  The  knowledge  that  is  gathered  can 
come  from  the  literature,  from  experts,  and  from  experiments  that  are  run  based  on  the  need  to 
acquire  specific  information  in  order  to  build  the  model.  Thus,  such  models  become  a  focal 
point  for  all  that  is  known  about  an  area,  and  an  exploratory  tool  for  gaining  further 
understanding  (Figure  1).  BioFusion  models  can  be  used  both  to  drive  in  vitro  and  in  vivo 
experimentation  as  well  as  to  serve  as  a  repository  for  the  results  that  are  acquired  from  such 
experiments.  In  this  way,  the  model  can  grow  with  the  knowledge  and  continue  to  provide  novel 
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insights  and  new  ways  of  looking  at  the  data,  as  well  as  help  prioritize  in  vitro  and  in  vivo 
studies. 


Literature 


i 


BioFusion® 
Model 


Experimental 
Data 


^  Expert  Opinion 
I  from  investigators 
I  in  the  group 


Figure  1  A  BioFusion®  model  incorporates  knowledge  from  literature,  experimental  data,  and  expert  opinion 

from  our  investigators 


The  development  of  a  BioFusion  model  involves  a  series  of  steps  centered  on  the 
acquisition  of  the  detailed  knowledge  required  to  build  a  model  that  exhibits  the  appropriate 
expected  behaviors.  The  BioFusion  model  development  process  imposes  a  methodology  on  how 
to  explore  a  given  area  targeted  for  modeling  and  how  to  analyze  and  synthesize  what  is  found 
into  a  coherent  picture,  or  theory,  of  the  process  being  studied.  A  key  attribute  of  a  BioFusion 
model  development  effort  is  a  series  of  meetings  where  the  information  that  has  been  collected  is 
sorted  through  and  examined  by  a  team  of  individuals  involved  with  the  development  of  the 
model.  This  team  of  individuals  must  include  the  experts  in  the  field  of  application  (such  as  the 
various  aspects  of  cellular  activity  in  the  eye  post  exposure  to  a  stimulus  such  as  a  laser),  the 
individuals  targeted  as  the  end  users  of  the  model,  and  the  computational  modeling  member(s) 
who  will  be  implementing  the  computer  model.  These  meetings  are  structured  and  managed  by 
the  Lead  Computational  Modeling  team  member  who  is  expert  in  this  knowledge  acquisition 
process  and  who  facilitates  and  directs  the  meetings  to  cover  and  analyze  what  is  needed  to  build 
the  model  in  the  selected  area  in  an  orderly  and  rigorous  manner. 

Developing  large,  complex  computer-based  simulations,  such  as  BioFusion  models,  is 
similar  to  developing  large,  complex,  knowledge -based  or  intelligent  systems.  In  both  cases, 
very  specific,  detailed  knowledge  must  be  acquired  and  represented  in  a  working  software 
program,  and  the  system  behavior  must  be  evaluated  and  accepted  by  an  expert  or  set  of  experts. 
In  both  cases,  the  development  of  a  working  software  system  is  more  a  problem  of  exploration 
than  of  implementation  (Fink  and  Herren,  1997),  because  system  requirements  are  not  easily, 
clearly,  or  completely  definable  at  the  start.  In  such  situations,  specification  and  implementation 
are  necessarily  intertwined  (Swartout  and  Balzer,  1982).  The  key  methodology  used  in  the 
development  of  knowledge -based  or  intelligent  systems  is  the  knowledge  engineering  process. 
The  Principle  Investigator  on  the  project.  Dr.  Fink,  has  had  extensive  experience  in  the 
knowledge  engineering  process,  and  utilizes  a  software  development  methodology  known  as  the 
Spiral  Development  Methodology  (Figure  5),  a  common  approach  in  software  engineering  of 
complex  computer  systems.  The  spiral  software  development  process  is  highly  iterative  and 
capitalizes  on  the  use  of  rapid  prototyping  to  define  the  behavior  of  the  final  system.  The  width 
of  the  spiral  at  each  point  represents  the  size  of  the  growing  software  system.  Thus,  the  approach 
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advocates  starting  out  small  and  adding  capabilities  incrementally  as  understanding  of  the 
problem  increases. 

Figure  2  demonstrates  this  interactive  and  iterative  process.  The  general  process  consists 
of  four  major  phases  of  development:  acquisition,  formalization,  implementation,  and 
verification/validation.  As  applied  to  simulation  models,  the  development  phases  are  knowledge 
acquisition,  knowledge  formalization,  model  implementation,  and  model  behavior 
verification/validation.  The  iterative  process  used  in  BioFusion  model  development  encompasses 
all  four  stages.  At  each  point  in  a  model  development  effort,  the  current  focus  defines  what 
knowledge  needs  to  be  collected,  formalized,  implemented,  and  validated.  The  initial  iteration 
generally  involves  only  knowledge  acquisition  and  formalization  and  produces  a  general 
description  of  the  system  and  listing  of  the  important  components  (e.g.  the 
cells/molecules/behaviors  involved  in  the  biological  process  of  interest).  This  information  is 
documented  in  a  set  of  knowledge  diagrams.  These  knowledge  diagrams  illustrate  graphically 
the  relationships  between  important  aspects  of  the  biology  being  modeled.  Depending  on  the 
targeted  level  of  the  model  to  be  developed,  the  relationships  may  be  between  cells,  between 
cells  and  various  cytokines/biologically  active  molecules,  or  between  various  chemicals  relevant 
to  the  key  internal  processes  to  the  cell.  In  addition,  all  references  used  are  documented  in  an 
electronic  library. 


Figure  2  The  spiral  software  development  methodology  used  to  develop  BioFusion®  models 


Subsequent  iterations  in  the  process  involve  defining  and  building  models  of  the 
components  and  then  validating  the  behavior  of  each  component  in  isolation.  Later  iterations 
involve  defining  and  implementing  interactions  between  components  and  validating  the  behavior 
of  groups  of  components  (e.g.  how  the  cells  interact  based  on  physiological  and  anatomical 
characteristics,  or  how  intra-cellular  cascades  interact  or  result  in  output  from  the  cell).  Final 
iterations  involve  global  impacts  on  the  system  components  and  validating  the  correctness  of  the 
global  changes  across  all  components  of  the  model.  In  each  iteration,  not  only  is  the  model 
expanded  and/or  refined,  but  the  knowledge  diagrams  and  the  reference  library  are  updated,  as 
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well,  to  reflect  the  new  knowledge/understanding.  In  this  way  a  BioFusion  model  is  fully 
documented  at  all  levels. 

The  BioFusion  approach  to  modeling  biological  systems  can  be  implemented  in  a  variety 
of  ways.  Usually,  however,  a  standard  commercial  off-the-shelf  (COTS)  model  development 
tool  is  used  to  generate  the  model.  The  COTS  tool  of  choice  has  been  a  general  modeling  tool 
called  Extend®,  by  Imagine  That!,  Inc.  Supporting  tools  used  for  the  knowledge  diagrams  and 
the  reference  library  are  Visio  and  EndNotes,  respectively.  These  tools  have  worked  very 
successfully  on  numerous  other  BioEusion  modeling  efforts  and  thus,  this  was  the  approach  used 
in  this  project,  as  well. 


2.0  SUMMARY  OF  WORK  PERFORMED 

A  BioEusion  model  of  the  effects  of  lasers  on  the  retina  was  developed  for  the  Air  Eorce 
Research  Labs’  Directed  Energy  Division.  The  work  focused  on  the  collection,  analysis, 
organization,  documentation,  and  implementation  of  information  particularly  on  the  retinal 
pigmented  epithelial  (RPE)  cell,  and  its  relationship  to  the  photoreceptor  cell. 

Over  the  24-month  course  of  the  project,  the  BioEusion  modeling  methodology  was 
iterated  through  a  number  of  times.  Ultimately  a  model  was  developed  that  focuses  primarily  on 
the  various  oxidative  and  photochemical  stress  mechanisms  that  take  place  inside  the  RPE  cell  as 
a  result  of  exposure  to  light.  The  model  is  primarily  one  of  intracellular  processes  to  the  RPE 
cell,  but  with  links  to  photoreceptor  influences,  as  well  as  other  environmental  effects. 

During  the  course  of  this  project,  a  set  of  knowledge  diagrams  covering  the  processes 
involved  with  photochemical  processes  and  oxidative  stress  in  the  RPE  cells  was  developed, 
along  with  an  electronic  reference  library.  In  total,  18  knowledge  diagrams  were  developed  and 
use  for  the  model,  along  with  two  knowledge  diagrams  related  to  details  on  the  visual  cycle  that 
were  only  used  as  a  higher  level  source  of  knowledge  for  the  model.  Hundreds  of  research 
articles  were  collected  in  support  of  the  knowledge  acquisition  process,  but  this  number  was 
narrowed  down  to  approximately  250  that  were  actually  used  in  the  model  development  and 
recorded  in  the  electronic  reference  library  (using  EndNotes).  The  rest  of  the  papers  collected 
and  analyzed  did  not  end  up  being  in  areas  that  became  the  focus  for  this  research  project, 
including  RPE  phagocytosis  of  photoreceptor  outer  segments,  the  interphotoreceptor  matrix,  and 
the  visual  cycle.  A  step  between  the  reading  of  the  reference  papers  and  the  development  of  the 
knowledge  diagrams  is  the  generation  of  relevant  “factoids”  about  the  various  topics  of  interest. 
These  are  generated  both  as  a  direct  result  of  reading  the  papers,  as  well  as  from  meeting  and 
discussing  the  information  in  the  papers.  These  factoids  are  recorded  in  a  set  of  Word 
documents  that  maintain  the  research  “history”  for  the  project.  Approximately  2,500  such 
factoids  were  recorded  in  the  process  of  developing  the  model  for  this  project,  covering  topics 
ranging  from  photoinitiation  and  redox  processes,  to  lipoprotein  damage  and  oxygen  control  in 
the  retina. 
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2.1  Overview  of  the  Laser-Retina  Model 

Figure  3  provides  a  top-level  knowledge  diagram  illustrating  the  areas  covered  in  the  the 
Laser-Retina  model.  This  top  level  knowledge  diagram  shows  the  relationships  between  the 
major  components  of  the  model.  Although  the  photoreceptor  and  visual  cycle  are  illustrated  in 
theis  diagram,  these  areas  were  not  implemented  in  the  actual  computer-based  model  in  detail. 
Figure  4  is  an  example  of  a  lower  level  knowledge  diagram,  one  at  the  actual  model 
implementation  level.  In  particular,  this  knowledge  diagram  illustrates  a  portion  of  the 
reduction- oxidation  processes  that  are  significant  to  the  RPE  cell’s  response  to  laser  exposure. 
The  model  contains  a  total  of  18  diagrams,  including  ones  for  photoinitiation,  the  melanin  effect, 
antioxidant  enzymes,  the  reduction  chain,  metabolism,  lipid  peroxidation,  DNA  damage,  and  the 
visual  cycle,  among  others.  The  emphasis  of  the  knowledge  content  is  on  oxidative  processes 
because  they  are  instrumental  to  the  photoreceptor  and  RPE  cells’  response  mechanisms  to  light 
exposure. 


Light-Retina  Biochemical  Interaction 
RPE  and  Photoreceptor 


Drafts: 

5/26/09 


Figure  3  The  top  level  knowledge  diagram  for  Laser-Retina  Model 
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Drafu: 


Figure  4  An  example  lower  level  knowledge  diagram  for  the  Laser-Retina  Model 


Figure  5  provides  an  example  of  what  the  Laser-Retina  Model  looks  like.  As  mentioned 
previously,  the  model  is  implemented  in  a  software  development  tool  called  ExtendSim,  by 
Imagine  That!,  Inc.  This  tool  is  a  visual  model  development  tool  where  the  modeling  blocks  can 
be  customized  to  meet  specific  modeling  needs.  Figure  5  illustrates  the  hierarchical  nature  of  the 
Laser-Retina  model,  which  is  supported  by  the  ExtendSim  model  development  tool.  The  figure 
shows  the  top  level  of  the  model,  along  with  the  next  lower  level  of  the  model  for  the  Redox 
Chain,  and  then  further  details  for  the  GSH/GSSG  portion  of  the  redox  chain. 

To  run  a  simulation,  attributes  of  a  light  exposure,  including  power,  spot  size, 
wavelength,  and  duration,  is  entered.  The  model  is  designed  initially  to  run  in  a  balanced  state 
with  nominal  visible  light  exposure.  Figures  6  and  7  show,  respectively,  a  portion  of  the  input 
screen  and  a  few  resulting  charts  from  a  simulation  run  where  not  only  ordinary  ambient,  visible 
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light  exposure  is  input,  but  also  a  laser  pulse  of  10  W,  with  laser  diameter  (spot  size)  of  1  cm, 
and  exposure  duration  of  10  seconds  is  introduced  at  timestep  500.  A  timestep  within  the  model 
is  set  to  be  equivalent  to  one  second  real  clock  time.  The  model  is  designed  to  run  for  the 
equivalent  of  48-72  hours  (2880  -  4310  minutes,  or  172,800-258,600  seconds),  but  it  can  be  set 
to  run  for  any  duration.  Actual  model  run  time  for  a  forty-eight  hour  simulation  takes  about  6-7 
seconds  on  a  standard  Windows  PC. 


1^'  ExtendSim 


Frie  Edit  Text  Library  Model  Database  Develop  Run  Window  Help 

□  if  #  ..«#*>•  100»  J  E'  -f  aa  ®  @  ^ '  [V  ^  ^  ii  ii  'S|  C  1  X:1296Y:58~|  0 


4' [1463]  Glutathione  <RPE2_DBJ.mox> 


o  'Laji.sj 


Figure  5  An  example  portion  of  the  Laser-Retina  Model  showing  the  hierarchical  nature  of  the  model 
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Figure  6  A  portion  of  the  input  screen  for  the  Laser-Retina  Model 
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Figure  7  Some  example  output  charts  of  the  results  of  running  the  Laser-Retina  Model 


The  effect  of  the  laser  pulse  input  can  be  seen  in  the  five  charts  shown  in  Figure  7.  The 
charts  (from  upper  left,  and  across  to  lower  right)  present  the  modeled  response  in  the  RPE  cell 
of  Vitamin  C  (asc/ascH),  some  reactive  oxygen  species  (including  superoxide,  hydrogen 
peroxide,  and  the  hydroxyl  radical),  NAD/NADH,  and  glutathione  (GSH/GSSG)  to  the  laser 
exposure.  A  portion  of  the  RPE  Cell  Health  chart,  which  provides  a  more  global  assessment  of 
the  RPE  cell’s  stress  level,  is  also  visible  on  the  right.  As  can  be  seen  in  these  charts,  the  model 
is  stable  prior  to  the  treatment,  has  a  response  for  approximately  two  to  two  and  half  hours  post¬ 
treatment,  and  then  re-stabilizes.  Thus,  this  level  of  treatment  produces  some  damage  to  the  RPE 
cell,  but  not  enough  to  actually  kill  it. 

The  goal  of  the  Easer-Retina  model  was  to  develop  a  tool  that  could  be  used  to  help 
explore  the  mechanisms  by  which  an  RPE  cell  becomes  damaged  and/or  dies  due  to  light 
exposure.  There  are  some  data  in  the  literature  regarding  laser-tissue  interactions  in  vitro  and 
even  though  there  appears  to  be  some  correlation  with  animal  studies  (Denton  et  ah,  2008, 
Schulmeister  et  ah,  2008),  direct  extrapolation  of  the  in  vitro  data  would  be  problematic  at 
present.  We  have  designed  the  model  to  allow  for  bifurcated  configurations  for  in  vitro  and  in 
vivo  environments.  In  the  retina,  RPE  cells  have  a  cuboidal  shape  with  the  melanosome  particles 
layered  just  beneath  the  microvilli  at  the  apical  surface.  The  in  vitro  configuration  excludes  the 
ocular  layers  anterior  to  the  RPE  layer  that  are  considered  as  optical  filters  in  the  in  vivo 
configuration,  and  it  accounts  for  the  flattened,  fried  egg-like  shape  of  RPE  cells  grown  on 
plastic  as  monolayers.  This  flattened  shape  exposes  cellular  organelles  such  as  the  nucleus  and 
the  mitochondria  to  greater  doses  of  incident  light  as  compared  to  the  RPE  cell  layer  in  vivo  for  a 
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given  laser  power  density.  Additionally,  the  in  vitro  scenario  does  not  account  for  the  complex 
relationships  between  retinal  RPE  cells  and  the  photoreceptor  cells  that  they  service,  such  as 
providing  needed  nutrients,  removing  waste  products,  recycling  retinal  via  the  visual  cycle,  and 
phagocytosing  shed  outer  segments.  It  is  possible  that  the  mechanisms  that  damage  and/or  kill 
an  RPE  cell  in  one  environment  are  quite  different  from  those  that  damage  and/or  kill  it  in  the 
other.  The  model  should  help  to  delineate  these  issues. 


2.2  Overview  of  the  Model  Database  and  Interface 

The  Easer  Retina  Model  contains  approximately  40  inputs  and  50  outputs,  representing 
the  various  molecules  or  ratios  that  are  tracked  by  time  step  throughout  the  course  of  a  model 
run.  As  a  result,  thousands  of  data  values  are  produced  and  recorded  during  each  simulation  run. 
These  data  are  initially  only  available  within  the  ExtendSim  modeling  and  simulation 
environment,  which  only  stores  the  results  from  a  total  of  four  simulation  runs  at  a  time  in  chart 
form  and  only  the  most  recent  run  on  table  form  with  actual  values.  Older  results  are  dumped  in 
chronological  order  as  new  simulation  runs  are  made.  In  order  to  make  data  available  for  analysis 
to  scientists  performing  research  beyond  these  few  simulation  runs,  the  data  must  be  exported 
from  the  Easer-Retina  Model  and  stored  into  a  database  for  later  access.  Thus,  a  database  and 
user  interface  was  developed  to  support  this  need.  This  “reporting”  system  offers  a  customizable 
interface  to  the  user  so  that  the  user  can  easily  select  the  data  they  need  to  create  a  report  in  the 
desired  report  format.  This  customizable  report  format  allows  for  the  selection  of  graphs/charts, 
if  wanted,  or  just  numerical  tables.  The  structure  of  the  report  system  includes  three  main  parts, 
1)  the  data  source  and  import  facility,  2)  a  customization  interface,  and  3)  the  reports. 

The  data  source  is  the  Easer-Retina  Model  itself.  Each  simulation  run  produces  a  pipe- 
delimited  (e.g.  “I”)  text  file  containing  the  value  of  all  inputs,  and  the  values  of  all  molecules 
available  in  the  model  on  each  time  step  as  output.  Each  file  generated  is  named  based  on  a  root 
name  and  a  sequence  number  so  that  the  results  of  each  simulation  run  is  stored  uniquely  and  not 
written  over.  Each  of  these  files  can  then  be  imported  into  the  model  reporting  tool  that  provides 
a  database  for  storage  and  an  interface  for  selecting  and  viewing  the  model  run  data.  The  model 
database  and  reporting  tool  is  implemented  in  Microsoft  Access.  Eigure  8  provides  an  example 
of  the  contents  of  the  exported  model  data  file;  Eigure  9  shows  the  first  screen  of  the  model 
database,  where  the  file  can  be  imported  into  the  database,  and  the  resulting  imported  data;  and 
Eigure  10  shows  the  user  can  select  the  particular  simulation  run  values  to  view.  Eigure  11 
shows  a  resulting  graph  that  plots  the  data  across  time  from  the  simulation  run.  This  chart  is 
similar  to  the  ones  produced  within  the  ExtendSim  model,  but  are  based  on  a  sampling  of  the 
data  every  minute  instead  of  every  second  to  save  storage  space. 
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outva2.txt  -  Notepad 

rsirwii;^ 

File  Edit  Format  View  Help 

0| 2. 1 11 1 0|0 10. 15 1  OHIO  10.45 110 110000 1104. 293  1 1000 10!0|0 1 808 1184 1492  I  500 115. 445 1 984. 555 1 200] 50 |100|10C  > 


13.1291371775518 

17.6771469812561 

17.6907138786033 

17.6848123277116 

17.6771004846198 

17.6724008441716 

17.6697320484568 


17. 668447180796 1  3. 82709116900191  HI 010116. 490021415682  I  5. 7818373481862  1 0149. 4700642470461 11099.  334763 


17.6681190945606 

17.6684670234842 

17.6693076607468 

17.6705223079106 

17.6720351122051 

17.6737986445002 

17.6757843654141 

17.6808094311847 

19.1710281705635 

17.7023693217441 

17.7500273085844 

17.7535614275678 

17.7570898378383 

17.7608597524186 

17.7648999943571 

17.7692075349147 

17.7737952650143 

17.7786753562228 

17.7838609874344 


17.7952058593338 

17.8013955190272 


4. 83568076198359  HI  010114. 8242494333151 1 01441. 392417529409 144 . 4727482999454 1 988. 283^ 
3.82676996277148111010115.1074339533182101214.905234895933145. 322301859954511007.1621=1 
3. 82619176775356  HI  010115. 1346648798272 14. 54244202133843 1 0|45. 4039946394817 11008. 97/1' ' 
3.  82631546117483  HI0I0H5.  376597505544  |4. 73573690480219 10 146. 1297925166321 11025. 106!  " 
3. 82668013198625  HI 0101 15. 6344803733966 14.97772363894001 1 0146.9034411201897  1 1042. 29{ 

3.  82690269487174  HI0I0H5.  90555875419251  5.  232312284 6423  10147.  716676262577611060.  370! 

3. 82702953432944  H 1 0 1 0 116. 1904992059187 | 5. 50012134135816  j  0 148. 57149761775611079. 366t 


3. 82710775237122 H 1010116. 804888710799816.07818905694287)0150.414666132399311120. 32! 
3. 82709253432971 H 1010117. 1358959831937 16. 3899357982771 |0| 51. 4076879495811 11142. 393C 
3.  8270541676742711)010117.483884933379516. 7178823301529810152.451654800138511165. 592 
3. 826998257459)1 10 10 117. 849737259650717. 06287238634945)01 53. 54921177895211189.982482 
3.82692838763823)1)0)0118.2343799561778)7.425793893866710)54.7031398685335)1215.6252 
3.  82684680141864)1)010118. 6387880354602)7. 80758173210283101 55. 9163641063806)1242. 58! 
3.82675485105812111010119.063987063247218.20922033140062)0157.1919611897417)1270.932 
3. 82652447320229 jl|0|0|19. 511728996408918.65744129717392)0158. 535186989226611300.781 
3. 71648369492631HI0|0|21. 0276202892732 19. 69441795586112 10163. 0828608678196)1401. 841 
3. 82561529782892jl)0|0|21. 706395431621)10. 8980743951022)0)65.11918629486311447. 09302 
3. 82331188133254 HI0|0|22. 3226558532413  111. 842850359232 ! 0 ) 66. 9679675597239 1 1488. 177C 
3. 82314165175178jl|0|0|22. 985602131746112.4824657407387)0)68. 9568063952379)1532. 3734 
3. 82297805051208HI0|0|23. 6841154996974113.156390371026210)71. 052346499092311578. 941 
3.  82280327558601H10|0|24.418945135525 113. 8652832309907 10)73. 2568354065748 11627. 929^ 
3.82261586438411110)0)25.1941957894559)14.6131862982092)0)75. 5825873683677 ) 1679. 613C 
3. 82241608897922)1)0 10126. 010890601196 115.401025213792810) 78. 0326718035881 11734. 059! 
3. 82220330849933)1)010) 26. 87119542943)16. 2308701188986 )0180.61358628829H791.413028C 
3.82197695901889)1)0)0)27.7773832523357)17.104894566298910183. 332149757007111851.82! 
3.82173642947544)11010128.7318384650223118.025377959041210186.195515395066911915.45! 


17. 78936620861 1  3. 82148106814357 11)0 10 [29. 7370615866113118. 9947099373969)0 [89. 2111847598338)1982. 47077 


3. 82121018648707 11)0)01 30. 7956740789093 120. 0153949064329 [0192.3870222367278  I  2053. 044 
3. 82092306151509 1 1 ) 0 | 0 ) 31. 9104232697408 [ 21.0900566591723 ) 0 | 95. 7312698092223 | 2127. 363 


17. 807951477282  I  3. 82061893718772  HI  0)0) 33. 0841873710142 122. 2214430846222 10) 99.  2525621130427  1 2205.  6124 
17.8148907167159)3.8202970252290911)0(0)34. 3199805796023 | 23. 4124309489788 | 0 1102. 959941738807 | 2287. 99? 
17. 8222309032393) 3. 81995650558898 11)0)0) 35. 6209582469391 (24. 6660307362581 [0)106. 862874740817 [2374. 73C 
17.8299903802491)3.81959652671704)1)010)36.9904221008353)25.9853915323113)01110.97126630250612466.02? 
17.838188164318513.81921620576412111010)38.4318255003612)27. 3738059336512 | 0 [115. 295476501083 | 2562. 123 
17. 8468439406336) 3. 81881462879281)1)0)0) 39. 9487787017429)28. 834714959749110)119. 846336105229)2663. 253 
17. 8559780568653  I  3. 81839085105881110)0)41. 5450541100329)30. 3717129444259)01124.635162330099)2769.6702 
17. 865611514433 1  3. 8179438974119)1 (0)0 [43. 2245914878414  I  31. 9885523786589)0)129. 673774463524  1 2881. 6394! 
17. 8757659562866) 3. 81747276285757)1)0)0)44. 991503088643  I  33. 6891486735354 10)134. 974 509265929)2999.433! 


Figure  8  Example  contents  of  the  text  file  output  from  the  Laser-Retina  Model 


Clr,  d 

Hon 


Database  Tools 


input_date  ♦ 
201008190001 
201008190002 
201008250001 
201008270001 
201008280001 
201008280002 
201008290001 
201009030001 
20101226^' 


nick_name 

bbbb 

eeee 

aaaa 

test8888&8 

TODAY 

ucuc 

tests 


02_used_metab  »  02_used_ph  »  photochem_  »  melan 


13.1291371775518 

17.6771469812561 

17.6907138786033 

17.6848123277116 

17.6771004846198 

17.6724008441716 

17.6697320484568 

17.6684471807% 

17.6681190945606 

17.6684670234842 

17.6693076607468 

17.6705223079106 

17.6720351122051 

17.6737986445002 

17.6757843654141 

17.6808094311847 

19.1710281705635 

17.7023693217441 

17.7500273085844 

17.7535614275678 

17.7570898378383 


2.1 

4.83568076198 

3.82676996277 

3.82619176775 

3.82631546117 

3.82668013199 

3.82690269487 

3.82702953433 

3.82709116900 

3.82710775237 

3.82709253433 

3.82705416767 

3.82699825746 

3.82692838764 

3.82684680142 

3.82675485106 

3.82652447320 

3.71648369493 

3.82561529783 

3.82331188133 

3.82314165175 

3.82297805051 


^  Main  Form 
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REPORT  SYSTEM 


Pate)  [^2_J  /  I  26  j  /  j  2010  | 

□  To  Replace  the  data  existing  in  the  database 


lnput.txt  I 


C:\DATA\AFRL_Modeling\Models\input2.txt 


\DATA\AFRL_Modelinq\MQdels\outva2.t 


Go  to  Report  Interface 


Browser  J 


Browser  J 
[  Import  ) 


0  17.4338849334 
0  17.8497372597 
0  18.2343799562 
0  18.6387880355 
0  19.0639870632 
0  19.5117289964 
0  21.0276202893 
0  21.7063954316 
0  22.3226558532 
0  22.9856021317 
0  73.6841154997 


6.71788233015298 

7.06287238634945 

7.4257938938667 

7.80758173210283 

8.20922033140062 

8.65744129717392 

9.69441795586112 

10.8980743951022 

11.842850359232 

12.4824657407387 

13.1563903710262 


Figure  9  Interface  for  importing  simulation  data  into  the  model  database  and  resulting  data  table 
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S)  i: 


input_date 

201008190001 

201008190002 

201008250001 

201008270001 

201008280001 

201008280002 

201008290001 

201009030001 

20101226^* 

'  3- 


nick_name 

bbbb 

eeee 

aaaa 


vitro  or  viv 


TODAY 

ucuc 

tests 


02_used_metab  -  02_used_ph  »  j 


13.1291371775518 

17.6771469812561 

17.6907138786033 

17.6848123277116 

17.6771004846198 

17.6724008441716 

17.6697320484568 

17.668447180796 

17.6681190945606 

17.6684670234842 

17.6693076607468 

17.6705223079106 

17.6720351122051 

17.6737986445002 

17.6757843654141 

17.6808094311847 

19.1710281705635 

17.7023693217441 

17.7500273085844 

17.7535614275678 

17.7S70R9837RW3 


2.1 

4.83568076198 

3.82676996277 

3.82619176775 

3.82631546117 

3.82668013199 

3.82690269487 

3.82702953433 

3.82709116900 

3.82710775237 

3.82709253433 

3.82705416767 

3.82699825746 

3.82692838764 

3.82684680142 

3.82675485106 

3.82652447320 

3.71648369493 

3.82561529783 

3.82331188133 

3.82314165175 

^■87797805051 
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Figure  10  Interface  for  selecting  data  from  the  model  database  to  present  in  a  chart  or  charts 
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Figure  11  An  example  chart  of  selected  data  from  a  Laser-Retina  Model  simulation  run 
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3.0  CONCLUSIONS 


The  BioFusion  modeling  approach:  1)  is  hierarchical,  2)  is  modular,  3)  can  model 
sufficient  depth,  4)  can  model  sufficient  breadth,  and  5)  is  synergistic.  This  technique  allows  for 
sufficient  depth  at  the  cellular  and  sub-cellular  levels  to  model  potential  interventions  at  their 
source,  while  at  the  same  time  covering  sufficient  breadth  of  the  relevant  biology  to  model 
events  at  the  organ,  system,  and  organism  levels.  Because  of  the  hierarchical  and  modular  nature 
of  the  resulting  models,  the  simulations  are  sufficiently  complex  to  allow  for  many  new,  unique, 
and  revolutionary  insights  into  the  biology  being  explored.  BioFusion  models  can  be  used  both 
to  drive  in  vitro  and  in  vivo  experimentation  as  well  as  to  serve  as  a  repository  of  results  found  in 
the  literature.  In  this  way,  the  model  can  grow  with  the  knowledge  base  and  continue  to  provide 
novel  insights  into  mechanisms  and  interactions  among  cells  and  systems. 

The  project  reported  on  here  has  addressed  a  key  component  of  the  retina,  namely  the 
photoinitiation  and  oxidative  stress  response  of  the  RPE  cell  to  laser  exposure,  as  the  starting 
point  for  further  understanding  of  how  light  affects  the  eye.  Use  of  the  modular,  hierarchical 
approach  to  developing  the  model  will  allow  work  to  focus  on  the  exploration  of  the  effects  on 
the  RPE  cell’s  viability  initially,  and  later  to  expand  to  include  details  of  other  cell  types  and 
processes  found  in  the  eye,  in  response  to  light  exposure. 
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